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Abstract 

We study the nonlinear dynamics of a deformed Deoxyribonucleic acid (DNA) 
molecular chain which is governed by a perturbed sine-Gordon equation coupled 
with a linear wave equation representing the lattice deformation. The DNA chain 
considered here is assumed to be deformed periodically which is the energetically 
favourable configuration, and the periodic deformation is due to the repulsive force 
between base pairs, stress in the helical backbones and due to the elastic strain 
force in both the strands. A multiple scale soliton perturbation analysis is carried 
out to solve the perturbed sine-Gordon equation and the resultant perturbed kink 
and antikink solitons represent open state configuration with small fluctuation. The 
perturbation due to periodic deformation of the lattice changes the velocity of the 
soliton. However, the width of the soliton remains unchanged. 
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1 Introduction 

Deoxyribonucleic acid (DNA) plays an important role in the conservation and 
transformation of genetic information in biological systems [1]. Opening of 
base pairs in DNA double helix is related to functions like transcription and 
replication. Base pair opening via nonlinear molecular excitations has been 
understood by several authors [2,3,4,5,6,7,8,9,10] by proposing different mod- 
els. Among them the models proposed and used by Yomosa [3,4], as well as 
by Takeno and Homma [5,6] were based on rotation of bases in a plane nor- 
mal to the helical axis of DNA, and the nonlinear molecular excitations were 
governed by kink-antikink solitons. Following Takeno and Homma, recently 
several authors [11,12,13,14,15,16,17,18] studied soliton-like molecular excita- 
tions in DNA by taking into account the rotation of bases. In all the above 
studies, both the strands of the DNA double helix were considered as rigid 
lattices. However, in nature the force between purine bases in consecutive base 
pairs is repulsive, and this force is resisted by stress in the helical backbones 
of DNA and also, the main-chain torsion angle indicates that there are elastic 
strain forces in both the strands [19,20]. The dynamics of this non-rigidity of 
the strands gives rise to phonons which also play an important role in energy 
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transfer in biological systems. In a different context, Davydov [21,22] proposed 
a model for energy transfer in alpha helix protein molecules and he found 
that the propagation of molecular vibrations induce longitudinal sound waves 
(phonons), which in turn provides a potential well that prevents vibrational 
dispersion, and this coupled excitation propagates as a soliton without loss of 
energy along hydrogen bonding spines of the alpha helical protein. Thus, the 
study of nonlinear molecular excitations in DNA double helical chains coupled 
with phonons in strands or in other words the influence of non-rigidity of the 
strands in molecular excitations has become an important task which requires 
a detailed investigation. In this direction, recently, Xiao and his co-workers 
[23] studied the influence of longitudinal vibration on the soliton excitations 
in DNA double helix by considering the dynamic plane base rotator model 
of Takeno-Homma [5,6], and by including the longitudinal vibration and its 
coupling with hydrogen bonds and stacking. It was shown that the dynamics 
in this case is governed by a perturbed sine-Gordon equation in the contin- 
uum limit, which upon solving using the method of successive approximation 
by iterations gives soliton under first order approximation, which shows that 
the effect of longitudinal vibration of the lattice on sohton is small. However, 
they failed to find the variation of the soliton parameters such as velocity and 
width explicitly during propagation under iterations. Therefore, in the present 
paper, we study the nonlinear molecular excitations in DNA double helix with 
non-rigid elastic strands, by solving the dynamical equation using direct soli- 
ton perturbation theory, which provides the variation of velocity and width of 
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the soliton under perturbation in explicit analytical terms. The paper is organ- 
ised as follows. In section 2, we consider the model Hamiltonian for the above 
DNA double helix and derive the dynamical equations. In the continuum limit, 
the dynamical equations reduce to a perturbed sine-Gordon equation coupled 
with a linear wave equation representing the longitudinal lattice vibration and 
this is treated in section 3. In section 4, a multiple scale soliton perturbation 
theory is developed to investigate the effect of lattice deformation on the open 
state configuration of DNA represented in terms of kink-antikink solitons of 
the sine-Gordon equation. The results are concluded in section 5. 

2 Model and dynamical equations 

We consider the B-form of a DNA double helix with flexible strands, and 
investigate the nonlinear molecular excitations by considering a plane-base 
rotator model. In Fig. (la) we have presented a sketch of the DNA double 
helix with z-axis parallel to the helical axis. In the figure, S and S' represent 
the two complementary strands in the DNA double helix, and each arrow in 
the figure represents the direction of the bases attached to the strand and the 
dots between arrows represent the net hydrogen bonding effect between the 
complementary bases. In Fig. 1(b), we present a horizontal projection of the 
n^^ base pair in the xy-plane in which Qn and Q'^ denote the tips of the n*'* 
bases, and P„ and represent the points where the n*'^ bases are attached 
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to the strands S and S' respectively. 



The conformation and stabihty of DNA double helical molecular chains are 
mainly determined by the stacking between the adjacent bases in the strands 
and the hydrogen bonds between the complementary bases. The Hamiltonian 
under plane base rotator model involving stacking and hydrogen bonds in 
terms of the rotational angles 0„ and 0^ of the n*'* bases (see Fig. 1(b)) in the 
case of rigid strands as proposed by Yomosa [3,4] and further developed by 
Takeno and Homma [5,6] is written as 



n 

-a[l - COS(0„ - (f)'J]] , (1) 



where / is the moment of inertia of the nucleotides around the axes at P„ and 
in the strands S and S' respectively and thus, the first two terms represent 
the kinetic energies of the rotational motion of the n*^ nucleotide bases. In 
Eq. (1) overdot represents time derivative. Further, the terms proportional to 
J in Hamiltonian (1) represent the stacking energy between the n*'* base and 
its nearest neighbours in the strands S and S' and a represents a measure 
of the interstrand interaction or hydrogen bonding energy between the com- 
plementary bases respectively. Pople's formula in which the mean energy of 
the distorted hydrogen bonds is approximately represented in the above form 
[24]. However, in nature DNA strands are not rigid but flexible and hence. 



we assume that the two strands deform elastically and the resultant phonons 
couple to the stacking and hydrogen bonds. Hence, the part of the Hamilto- 
nian corresponding to the phonon energy, and the energy due to its coupling 
with the stacking and hydrogen bonds is written as 



where p„ = My„ and = My'^. In Eqs. (2), y„ and y'^ represent the longitu- 
dinal displacements of the n*'* nucleotides from the equilibrium position in the 
two strands, and M is the uniform mass of the nucleotide. K is the longitudi- 
nal elastic constant along the double hehcal main chain. The stacking energy 
depends on the distance between the n^^ and (n -|- 1)*'^ base, and the strength 
of the hydrogen bonds depends symmetrically on the distance between the 
(n — and {n -\- 1)*'* bases. Thus, (3 and 7 measure the couphng strengths 
between phonon and stacking as well as hydrogen bonds respectively. The in- 
teraction Hamiltonian Hr-ph in Eq.(2b) is chosen to represent the change in 
stacking energy and hydrogen bonds energy caused by the change in the dis- 
placement of the nucleotides along the two strands. As we are going to study 
the dynamics in the low temperature and long wavelength limit, it is appro- 
priate to consider linear coupling of phonon to the stacking and hydrogen 
bonds. Now, using the Hamiltonians (1), (2a) and (2b), the total Hamiltonian 




(2a) 




J2 [P{{yn+i - yn)[l - cos((/)„+i - 0„)] + (y^+i - y'„) 



n 



x[l - cos(0;+i - - 7(2/^+1 - 2/n-i)[l - cos(0n - 0;)] (2b) 
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H for the system is written a,s H — Hr + Hph + Hr-ph and the corresponding 
Hamilton's equations of motion take the form 

lipn =[J + PiVn+l - Un)] sin(0„+i - 0„) - [J + /3(y„ - l/n-l)] 

X sin(0„ - + [a + 7(i/„+i - Vn-i)] sin(0„ - 0^), (3a) 

I(P'n =[J + PiVn+l - Vn)] sin(0^+i - (p'J - [J + (3{yn - Vn-l)] 

X sin((/>^ - + [a + 7(l/n+i - 1/n-i)] sin(0; - 0„), (3b) 

My„ = 2K{yn+i - 2y„ + y„_i) - /?[cos(0n+i - 0n) - cos(0„ - 0„-i)] 

+7[cos(0n+i - (f)'^^^) - COs((/)n_i - (/)^_i)], (3c) 

My; = 2X(y;+i - 2y'^ + y;_i) - /3[cos(0;+i - 0^) - cos(0; - ct>'^_^)\ 

+7[cos((/)„+i - 0;+i) - cos((/)„_i - (/);_i)]. (3d) 

Eqs. (3a-3d) describe the dynamics of DNA with deformable strands at the 
discrete level by considering the dominant angular rotation of bases in a plane 
normal to the helical axis of the DNA, and ignoring all other small motions 
of the bases combined with longitudinal motion of the nucleotides. 

3 Soliton and beise pair opening 

It is expected that the difference in the angular rotation of bases with respect 
to neighbouring bases along the two strands in DNA namely (0„±i — 0n) and 
(^n±i ~ ^n) small [5,6]. Also, as the length of the DNA chain is very 
large due to the presence of large number of bases compared to the distance 
between the neighbouring base pairs, we make a continuum approximation 
by introducing two fields of rotational angles 0n(^) — 0(2^, ^), 4''n{'^) ~^ 4>'{^j ^) 
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and two fields of longitudinal displacement — > y{z, t) and — > y'{z, t) 
along the strands where z = na with /, the lattice parameter. Also, we make 
the expansions for = (l)(z,t) ± + ± + and similar 

expansions for (p'n±i, yn±i and y'„-i-i- Thus, in the continuum limit under small 
angular rotation of bases Eqs. (3) upto O(a^) become 



1 7 

(t>ii^(t>zz - 2 sin((/> - 00 +e[/3(yA)z + 22/^si"(<^ - 0')], (4a) 

0;-,- = C - ^ sin(0' - 0) + e[/3(|/>;), + ||/; sin(0' - 0)], (4b) 

yti = v^yzz, (4c) 

y££ = ^'z/L, (4d) 



where e = j, = and the suffices t and z in Eqs. (4) represent partial 
time and spatial derivatives and the rescaled a is dimensionless. While writing 
the above equations we have chosen a — —\Ja? and also, the parameter 7 is 

2 

rescaled as 7 = Further, before writing Eqs. (4a-d), we have divided the 
full equations by J a? and rescaled the time variable as f = ^t. It is more 
convenient to describe the transverse motion of the bases in DNA strands 
in terms of the center of mass co-ordinates. For this, we rewrite Eqs. (4) by 
subtracting and adding the first two and the last two equations respectively. 
Further, to commence the open state configuration of DNA, the two com- 
plementary bases are expected to rotate in opposite directions and both the 
strands are assumed to vibrate in the same direction so that 4>' = —cj) and 
y' — y. Under these conditions, we obtain 



- + sin * = e[P{y,^,), + ^y, sin , 



(5a) 
(5b) 



where ^ = 20. Eqs. (5) describe the dynamics of bases under a plane-base 
rotator model of DNA double helical chain with the deformed strands. The 
terms proportional to /3 and 7 in the right hand side of Eq. (5a) represent the 
coupling of phonon to the stacking and hydrogen bonds respectively. 
When e = 0, Eqs. (5a) and (5b) are decoupled, and Eq. (5a) reduces to the 
completely integrable sine- Gordon equation which admits kink and antikink- 
type of soliton solutions, and hence we call Eq. (5a) in its present form as a 
perturbed sine-Gordon equation. The integrable sine-Gordon equation (e = 0) 
was originally solved for N-soliton solutions using the most celebrated Inverse 
Scattering Transform (1ST) method by Ablowitz and his co-workers [25] . The 
kink and antikink one soliton solution of the integrable sine- Gordon equation 
(Eq.(5a) when e = 0) can be written as 

'if{z,t} — 4:a,Tctanexp[±m{z — vi}], — Vl — v'^. (6) 

In Eq.(6), while the upper sign corresponds to kink soliton, the lower sign 
represents the antikink soliton. Here, v and are real parameters that 
determine the velocity and width of the soliton respectively. The kink and 
antikink one soliton solutions as given above are depicted in Figs. 2(a) and 
2(b). The kink-antikink soliton of the sine-Gordon equation describes an open 
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state in DNA double helix which is schematically represented in Fig. 2(c). In 
this figure the base pairs are found to open locally in the form of kink-antikink 
shape in each strand and the opening is found to propagate along the direction 
of the hehcal axis. 

Eq.(5b) is the well known one-dimensional linear wave equation which admits 
wave solution in the form y — f{z — vi) + g{z + vt), where / and g are arbitrary 
functions. Now, the problem boils down to solving the perturbed sine-Gordon 
equation (5a) after using the wave solution 'y' obtained by solving Eq. (5b). 

4 Effect of elastic deformation of strands on base pair opening 

4-1 A perturbation approach 

When the phonon due to elastic deformation of the strands is coupled to the 
DNA molecular excitations, it is expected to perturb the kink and antikink 
solitons in DNA which correspond to the open state configuration. It is further 
expected that the perturbation due to phonon coupling modifies the shape, 
width and velocity of the soliton as it propagates along the helical chain. In or- 
der to understand this, we solve Eq. (5a) using a suitable perturbation method. 
One of the most powerful techniques in dealing with perturbed soliton is the 
soliton perturbation theory which is based on the 1ST method. However, as 
the method is very sophisticated it is very difficult to use the same in several 
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cases. In view of this, many authors used different types of direct methods to 
study sohton perturbation (see for e.g. refs. [26,27,28,29,30]). In the present 
paper, we use one such direct perturbation method to solve the perturbed 
sine-Gordon equation (5a) to understand the effect of phonon interaction on 
the open state configuration of DNA, which is also dealt in reference [30] in a 
different context, and also by the present authors recently while studying the 
nonlinear molecular excitations in an inhomogeneous DNA [17]. The procedure 
we adapt here is based on the derivative expansion method to linearize the 
perturbed sine-Gordon equation in the coordinate frame attached to the mov- 
ing frame. The parameters of the kink-antikink soliton are assumed to depend 
on a slow time scale in order to eliminate the secular terms. The linearized 
equations will be solved using the method of separation of variables which will 
be ultimately related to a generalized eigenvalue problem, the eigenfunctions 
of which form the bases of the perturbed solution. In the following we use the 
above approach to find the perturbed soliton solution of Eq. (5a). 

4-2 Linearization of the perturbed sine-Gordon equation 

In order to study the effect of perturbation due to phonon interaction on the 
soliton, the time variable i is transformed into several variables as tn — eH, 
where n=0, 1, 2,... and e is a very small parameter. In view of this, the time 
derivative and ^ in Eq. (5a) are replaced by the expansions ^ = ^ + ^ ^ + 
+ ■■■ * = e^'(^) e^*(^) ... and we equate the coefficients of 
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different powers of e. Thus at 0{€^^^) we obtain 



*So-*l°;+sm*(°)=0, (7) 

for which the one sohton solution takes the form {z,to) — 4arc tan exp C, C = 
±mo(^ — 0, Cto — vo, where vq is the velocity of the soliton in the to time scale. 
Due to perturbation, the soliton parameters namely m and ^ are now treated 
as functions of the slow time variables to,ti,t2, ■■■■ However, m is treated as 
independent of to. The equation at 0{€^^^) is of the form 

^rc - *CC + (1 - 2sech2C)^(i) = F(^)(C, r), (8) 



where 



F(^) = 2/3 [y^sechC]^ + 2 b^y^ tanh (sechC 

+4^0 + (m^^tj — J tanh C sech^. (9) 

While writing the above equation we have replaced sin by 26 tanh ^ sech^, 
where b — ±1, which can be derived using the unperturbed solution given 
below Eq. (7), and we have also used the transformation ( — m{z — vto) and 
to = to to represent everything in a co-ordinate system that is moving with 
the soliton. Further, we have used another set of transformations given by 



T — — -^i^^ and C = ^ for our later 



convenience. 



The solution of Eq. (8) is searched by assuming ^'(^^(C,t) = X{()T{t) and 
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F(i)(C,t) = X^{C)H{t). Substituting the above in Eq. (8), we obtain 



^[Xcc + (2sech\ - 1)X] = 1[T, - H{t)] 



(10) 



C 



In Eq. (10), while the left hand side is independent of r, the right hand side 
is independent of the variable (. Therefore, we can equate the left and right 
hand sides of Eq. (10) to a constant, say Aq and write 



Thus, the problem of constructing the perturbed soliton at this moment turns 
out to be solving Eq. (11) by constructing the eigenfunctions and finding the 
eigenvalues. The first part of Eq. (11) is a generalized eigenvalue problem, 
which is not a self-adjoint eigenvalue problem and differs from the normal 
eigenvalue problem, with in the right hand side instead of X. For solving 
the eigenvalue problem, we consider it in a more general form by writing 



where A is the eigenvalue. In order to solve Eq.(12) for X, we also consider 
the following eigenvalue problem. 



X^^ + (2sech^C - 1)^ = AqX^, - AqT = H{t) . 



(11) 



LiX = XX, Li = % + 2sech2C - 1, 



(12) 



L2X = XX, 



(13) 
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where L2 is to be determined. Now, by combining the above two eigenvalue 
problems we get L2L1X = X^X, L1L2X = X^X. From these expressions we 
conclude that L1L2 is the adjoint of L2L1 and also X and X are expected to 
be adjoint eigenf unctions. Hence, by solving the coupled eigenvalue problem 
we can find the eigenfunction X. Here Li is known and is given in Eq. (12), 
but the operator L2 is still unknown. So, by experience we choose L2 — + 
6sech^^ — 1. 

Now, in order to find the eigenfunctions by solving Eqs. (12) and (13) we 
choose the eigenfunctions as 



X(C,A;)=p(C,A;)e^^ 1(C, A;) = ?(C, ^)e^^ (14) 

where k is the propagation constant. On substituting the above in Eqs. (12) 
and (13) in the asymptotic limit, we obtain the eigenvalue as A = —(1 + k"^). 
In order to find the eigenfunctions, we expand p(C, k) and q{(, k) as 



/> 7 \ sinhC 1 sinhC 1 x 

coshC cosh C cosh C cosh C 



i\ A_ sinhC ^ 1 ^ sinhC ^ 1 ^ 



where pj and g^, j=0,l,2,... are functions of k to be determined. Substituting 
Eqs. (14), (15a) and (15b) in Eqs. (12) and (13) and collecting the coefficients 

'-'^ 1' coshC ' cosh^c '"" ^ simultaneous equations. On solving those 

equations by assuming pj = qj — for j > 3, we obtain the eigenfunctions as 
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(1 - fc2 _ 2iA;tanhC) 

v^(l + P) ^ ' 
(1 - - 2iA;tanhC - 2sech\) 



(16b) 



(16a) 



On comparing Eqs. (16a) and (16b) we can write X((, k) = Now, using 



this in the right hand side of Eq. (12) and comparing the resultant equation 
with Eq. (11), we obtain Aq = ^'•^"^^ 

The second part of Eq. (11) is a hnear inhomogeneous differential equation 
and it can be solved using known procedures [17]. The solution reads 



The first order correction to the soliton can be computed using the following 
expression. 



Here X{(,k) and T{t, k) are known continuous eigenfunctions which are given 
in Eqs. (16a) and (17). However, the discrete eigenstates Xo,Xi and To,Ti 
are unknown. Xq and Xi are the two discrete eigenstates for the discrete 
eigenvalue A = and these states can be found out using the completeness of 
the continuous eigenfunctions as 



*(^)(C, r) = / X{C, k)T{r, k)dk + ^ X^{QT^{t) 
■1-°° i=0,l 



(18) 



Xo(C)=sechC, Xi(C)=CsechC. 



(19) 
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In order to find Tq and Ti, we substitute Eq.(18) in Eq.(8) and multiply by 

Xo(C) and ^i(C) separately, and after using the orthonormal relations, we get 



/oo 
FW(C,r)Xo(C)ciC, (20a) 
-oo 

/oo 
F«(C,r)Xi(C)ciC. (20b) 
-oo 

As F^^\(, t) given in Eq. (9) does not contain time r explicitly, the right hand 
side of Eqs. (20a) and (20b) are also independent of time, and hence they give 
rise to secularities and the nonsecular conditions can be written as 

/oo 
F^'\C,r)Xo{C)dC = 0, (21a) 
-oo 

/oo 
FW(C,T)Xi(C)dC = 0. (21b) 
-oo 

On substituting the above expressions in Eqs. (20a) and (20b), we choose 
Ti(r) = and obtain Tq{t) = C, where C is a constant, which has to be 
determined. For this, we integrate Eq. (20b) and obtain 

ro(r) = fl4^ r dC CF^'\Ct)X,{0. (22) 

Z J—oo 



4-3 Variation of soliton parameters 



In order to find the first order correction, we need to evaluate the eigenstates 
explicitly for which we need the values of m^^ and which can be found from 
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the nonsecularity conditions, by substituting the values of F^^\(, t),Xo{() and 
Xi{Q respectively from Eqs. (9) and (19). The results give the time evolution 
of the inverse of the width (m) and the velocity {^ti ) of the soliton as 



1 f°° 

mt^^-— {(3 [yi; sechQ + b^y^ tanh (sechC) sechC^C, (23a) 

ZVq J —oo 
1 f°° 

Cti = ~ 2mHo J- ^^c^Clc + ^^^C CsechOCsechC^C. (23b) 

In order to evaluate the integrals found in Eqs. (23a) and (23b) explicitly, 
we have to substitute the value of 'y' which we have found by solving Eq. 
(5b). We consider the most general and meaningful wave solution of Eq. (5b) 
suitable for the problem as the periodic function |/ = sin^. At this point, it is 
worth mentioning that Dandoloff and Saxena [31] realized that in the case of 
an XY-coupled spin chain model which is identifiable with our DNA double 
helical chain model, the ansatz cn{(, k) with the limit k — > 0, energetically 
favours the periodically deforming spin chain. Hence by substituting = cos C 
in Eqs. (23a) and (23b) and on evaluating the integrals we obtain 



7r[7r/3-b7(4-7r)] 



The parameters m and ^ can be written in terms of the original variable t as 



eirlPir — 67(4 — tt)] , , 

™ = m„, (i^v = v„+ » . (25) 
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where 1/mo is the initial width of the sohton and vq is the uniform velocity of 
the soliton in the unperturbed limit. The first of Eq.(25) says that, the width 
(m~^) of the soliton remains constant. However from the second of Eq. (25), 
we find that the velocity of the soliton gets a correction. It is observed that the 
correction in velocity depends on the nature of (3 and 7 which can be either 
positive or negative for 6 = ±1. First, we consider the case corresponding to 
6 = 1. In this case, when [/^tt — 7(4 — tt)] > 0, the velocity of the soliton gets 
a positive correction and hence soliton may propagate along the DNA chain 
without forming a bound state. On the other hand, when [/?7r — 7(4 — tt)] < 0, 
the phonon due to lattice deformation either slows down the soliton or the ve- 
locity of the soliton remains unaltered. Finally, if the initial uniform velocity 
of the soliton before switching on the perturbation due to elastic deformation 
takes the value = ^^^^(4^~)l\g , the soliton is stopped by the deformation. 
The stability of the soliton is guaranteed in all the above cases. A similar 
argument can be made in the case of 6 = — 1 with [Pn — 7(4 — tt)] replaced by 
[/37r + 7(4 — 77)]. Recently Yakushevich et al [14] and Salerno [11] investigated 
the interaction of soliton with periodic sequence (periodic inhomogeneity) , 
and the results have very close analogy with our results here. It was shown by 
them that soliton can easily propagate along DNA without forming a bound 
state. It may also be noted that, Zhang et al [32] obtained similar results in the 
case of resonant kink impurity interaction and kink scattering in a perturbed 
sine-Gordon model. In a recent paper, Hwa et al [33], while studying the ther- 
modynamic and dynamic behaviours of twist induced denaturation bubbles in 
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a long, stretched random sequence of DNA using statistical mechanical mod- 
els, has shown the localization and delocalization of bubbles along the DNA 
chain. Finally, Eq. (25) is also similar in form to our recent results of pertur- 
bative analysis in the case of an inhomogeneous DNA [17,18]. Thus, we can 
say that the lattice deformation gives rise to inhomogeneity in the DNA chain. 



4-4 First order perturbed soliton 



Now, we exphcitly construct the first order correction to the one soliton 
by substituting the values of X{(,k), Xo{(), Xi{() and T{T,k), Tq{t) from 
Eqs. (16a), (19) and (17), (22) and that oi F^^\C,t) from Eq. (9) and use the 
values of rrit-^ and from Eqs. (24) in Eq. (18) to get 



*«(C,to) = - 

TT 



dk 



;i -A;2 -2iA;tanhC)e**^^ 



oo (1 + F)3' 

/oo 
dC'{l -k'^ + 2ik tanh C') [P sin C + {{P - &7) cos C' 
-oo 



TT 



7r/3 - 67(4 - tt))} tanhC']sechC'e 



f-ikC 



X 



+ (1 + T;)sechC 



X 



/oo ^ 
dC C '[/SsinC' + {(/?- 67) cosC' - « (tt/? - 67(4 - tt))} 
-00 o 



X tanh C'Jsech^C 



(26) 



We evaluate the integrals in Eq. (26) by finding the values of the residues at 
poles of different orders using residue theorem (for details see [17,30,34]). After 
lengthy algebra and some approximations the explicit form of the perturbed 
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kink (upper sign)-antikink (lower sign) one soliton solution in terms of the 
original variables is written as 



"^{Zjto) ~4 arc tan exp [mo (-2 — t'o^o)] 
e7i[n(3 — 67(4 — tt)] 



m{v-l) + 2vto sech[m{z - vto)]. (27) 



Having found ^{z,to) we find (f){z,to) using the relation = f- and plot the 
same in Figs. 3(a,b) by choosing (3 = j = l,b — —1 and vq — 0.4. Prom 
the figures, we observe that the lattice deformation introduces only small 
fluctuations in the form of periodic oscillations closely resembling the shape 
of the lattice deformation in the width of the soliton (see Figs. 3(a,b)). We 
have schematically represented this in Fig. 3(c), where the dotted line along 
the strands (lattice) represent the periodic deformation of the lattice. It shows 
that the lattice deformation in DNA does not affect opening of bases in DNA 
double helix. 



5 Conclusion 

In this paper, we studied the effect of phonon interaction on base pair opening 
in DNA by considering the dynamic plane base rotator model. The dynamics 
of this model in the continuum limit gives rise to a perturbed sine-Gordon 
equation coupled with a linear wave equation representing longitudinal lattice 
vibration, which were derived from the Hamiltonian consisting of the stacking 
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energy, hydrogen bonding energy, energy corresponding to the lattice defor- 
mation and its couphng with the stacking and hydrogen bonding energy. In 
the unperturbed hmit, the dynamics is governed by the kink-antikink sohton 
of the integrable sine-Gordon equation which represents the opening of base 
pairs in DNA without lattice deformation. In order to understand the effect 
of lattice deformation on the base pair opening, we carried out a perturbation 
analysis using multiple-scale soliton perturbation theory. From the results of 
variation of soliton parameters we observe that when the DNA lattice de- 
forms in a periodic way, the width of the soliton remains constant. However, 
the velocity of the soliton increases or decreases or remains uniform or even 
the soliton stops depending on the values of the coupling strengths (3 and 7. 
Interestingly, the soliton in all the above cases are found to be stable. From 
the results of the perturbed sohton we observe that the periodic lattice defor- 
mation introduces fluctuation in the width of the soliton. However, there is 
no change in the topological character of the soliton in the asymptotic region. 
The above dynamical behaviour may act as energetic activators of the enzyme 
transport during the process of transcription in DNA. 
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Fig. 1. (a) A schematic representation of the structure of B-form DNA double 
hehx. (b) A horizontal projection of the n^^ base pair in the xy-plane. 
Fig. 2. (a) Kink and (b) antikink one soliton solutions of the sine-Gordon 
equation (Eq. (5a) when e = 0). (c) A sketch of the formation of open state 
configuration in terms of kink-antikink solitons in a DNA double helical chain. 
Fig. 3. The perturbed (a) kink-sohton and (b) antikink-sohton with (3 — ^ — 
1.0,6 = —1.0 and Vq = 0.4. (c) A sketch of the open state configuration in 
DNA with small fluctuations and periodic deformation in the lattice. 
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Fig. 2. (a) Kink and (b) antikink one soliton solutions of the sine-Gordon equation 
(Eq. (5a) when e = 0). (c) A sketch of the formation of open state configuration in 
terms of kink-antikink sohtons in a DNA double helical chain. 
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(C) 

Fig. 3. The perturbed (a) kink-soliton and (b) antikink-soliton with 
/? = 7 = 1.0,6 = —1.0 and vo = 0.4. (c) A sketch of the open state configura- 
tion in DNA with small fluctuations and periodic deformation in the lattice 
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